Lattice animals and the Percolation model under rotational 

constraint 



Indrani Bose 

Department of Physics 
Bose Institute 
93/1, A. P. C. Road 
Calcutta- 700 009, India. 

Abstract 

The effect of rotational constraint on the properties of lattice mod- 
els like the self-avoiding walk, lattice animals and percolation is dis- 
cussed. The results obtained so far, using a variety of exact and 
approximate techniques, are described. Examples of the rotational 
constraint in real systems are also given. 

I. Introduction 

Lattice models describe a wide variety of physical, chemical and biological 
systems. The systems either have an underlying lattice structure or are stud- 
ied on a lattice for convenience. In the latter case, only those properties of 
the system are studied which are independent of the microscopic structural 
details of the system. In this article, we will discuss three lattice models: 
the self-avoiding walk (SAW), lattice animal (LA), and percolation. We will 
mainly discuss the effect of rotational constraint on the various configura- 
tional and connectivity properties of the models. 

A SAW is a walk on a lattice with the constraint that a site once visited 
cannot be visited again. The SAW is a model of linear polymers in dilute 
solution [l| . A linear polymer can be pictured as a string of smaller chemical 
units called monomers. There is a one-to-one correspondence between each 
possible configuration of the polymer and a particular SAW on the lattice. 
The self-avoidance criterion is due to the excluded-volume effect of polymers 
which implies that no two monomers can occupy the same physical space. 

LAs are clusters of connected sites (site LA) or bonds (bond LA) embed- 
ded in a lattice. The LAs represent random clusters which can occur in a 
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variety of systems. In terms of polymers, one example is branched polymers 
in dilute solution. The percolation model is a model of disorder and has, in 
general, two versions, site and bond. In the site (bond) percolation model, 
the probability that a site (bond) of the lattice is occupied is p. For p less 
than a critical value p c , called the percolation threshold, no extended network 
of sites (bonds) spans the lattice. For p i p c , such a network exists with prob- 
ability one. The percolation threshold p c has a unique value for an infinitely 
large lattice. The value of p c is different for different lattices. The percolation 
transition at p = p c is a phase transition, often called a geometrical phase 
transition. In one phase ( p I p c ) , an infinite cluster exists and in the other 
phase ( p j p c ), there is no infinite cluster. In the vicinity of p = p c , the 
percolation model exhibits critical phenomena analogous to those exhibited 
by a second-order thermodynamic phase transition. The percolation model 
has extensive applications in physics, chemistry and biology |2], |3|, |], |5j . 

We now show the connection between LAs and the percolation problem. 
We consider the site cluster problem. The generating function desfribing an 
ensemble of clusters of all sizes is written as 

G{x,y) = Osj^y* (1) 

s,t 

g s j is the number of distinct LAs of size s and t perimeter sites, each animal 
being weighted by a factor x s y l . The size s of the animal is equal to the 
number of sites in the cluster. In the percolation problem, each cluster site 
is occupied with probability p and a perimeter site being an unoccupied site 
has the probability 1-p of being so. Thus x = p and y = 1-p in Eq.(l).The 
corresponding generating function is 

G{p,q) = E# s ,tpV 

s,t 

= p a (2) 

s 

where q = 1-p and D s (q) = J2t 0s,t9*-The D s (q) are known as perimeter 
polynomials. For LAs, clusters of the same size have equal weight irrespective 
of the number of perimeter sites, so y = 1 in Eq.(l) and the generating 
function becomes 

G(x,l) = Y,9sX s (3) 
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where g s = J2t 9s,t is the total number of clusters of size s. Eq. (2) shows that 
percolation and LA models are related through the perimeter polynomials 
construction of which requires knowledge of LA numbers. Both generating 
functions (2) and (3) become singular as p — > p c and x — > x c , respectively, 
hence p c and x c are identified as critical points of the respective ensembles of 
clusters. Below the critical point small clusters are predominant in number, 
above the critical point the reverse situation is true. At the critical point, 
clusters of all sizes exist, a familiar characteristic of critical phenomena. In 
fact, the transition is a second-order transition and the generating function 
in each case plays the role of free energy in thermal critical phenomena. 
Analogous to thermodynamic quantities, cluster-related quantities can be 
defined for the geometric models. These show singular behaviour as the 
critical point is being reached. The singularity is of the power-law type with 
an associated critical exponent. For the SAW model, a generating function 
similar to Eq.(3) can be defined with s being the number of steps in the walk. 
For generating functions of the type (3), it can be shown || that every critical 
exponent characterising a singularity as x — > x c , say, is uniquely related to 
a critical exponent along the path s — > oo, i.e., s ~ | x — x c | -1 . In the 
asymptotic limit of s going to infinity, various walk or cluster properties scale 
with size. We will now describe two such relations for site LAs. Equivalent 
relations can be written down for SAWs. The first relation shows that the 
total animal (cluster) number g s scales as 

g s ~™ X s s- e (4) 

where A is a constant for a particular lattice and is known as the 'growth 
parameter' because asymptotically g s / g s -x —* A. The parameter 6 (critical 
exponent) is 'universal' in the sense that it has the same value for all lattices 
in dimension d, the value changes only if d changes. The second quantity 
we are interested in is the average radius of gyration defined as R s = (< 
J2t=i r i I s >) 1//2 where is the distance of a cluster site i from the centre of 
mass of the cluster and denotes average over all animals. As s — > oo, 
R s scales with s as 

R s ~ s v (5) 

where v is another critical exponent. In connection with relations like (4) and 
(5), it has been shown || that LAs and large clusters below the percolation 
threshold have identical scaling behaviour. For the SAW model, g s is the 
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number of distinct SAWs with s steps and R s is the root-mean-square end- 
to-end displacement for walks with s steps. 

The scaling relations (4) and (5) have been well-studied for both undi- 
rected and directed models. In the case of undirected SAW and LA, there 
is no constraint on connectivity. In the case of directed SAW and LA, there 
is a global directionality constraint. To give an example, starting from the 
origin and while occupying new sites, the animal can grow only upwards and 
towards the right. The undirected and directed lattice models are known 
to belong to different universality classes, i.e., the exponents 9 and v have 
different values. A new type of constraint ]7|, || [| [10] that has been proposed 
is the rotational constraint. In the presence of this constraint, a SAW or a 
LA can grow either in the same direction as it has been growing in or in a 
specific rotational direction, say, clockwise. Similar constraint can also be 
defined for the percolation model. In Section II, the SAW model under rota- 
tional constraint, termed the spiral SAW (SSAW), is described. In Section 
III, spiral LAs are considered. A brief description is also given of directed 
compact site LAs for which the generating function can be determined in 
an exact manner. In Section IV, the spiral percolation model is discussed. 
Section V contains some concluding remarks. 



II. Spiral SAW 

The model was first considered by Privman|7]]. Let us consider the SSAW 
problem on a square lattice. Each step of the walk may either point in the 
same direction as the preceding one or in a direction rotated by + tt/2 w.r.t. 
it. The global behaviour arising from the two constraints of self- avoidance 
and the bond angle restrictions is quite novel. In general, a SSAW consists of 
an outward spiralling part and an inward spiralling part (Fig.l). Each SAW 
is composed of horizontal and vertical segments. Let s n be the number of 
n-step SSAWs. In the asymptotic limit, exact, analytical expressions for s n 
have been determined [II, [12], [T^, O, |TJJ] in the cases of both the square and 



triangular lattices. We describe briefly the results obtained in Refs. [12-13] 
for the square lattice. 

Let the first step of the walk be vertically upwards. Each SSAW can 
be decomposed into a concatenation of two walks of type I and C by either 
a vertical line (Fig.2(a)), a horizontal line (Fig. 2(b)) or in some cases both 
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(Fig. 2(c)). I n denotes the set of n-step SSAWs whose last and third last seg- 
ments are of equal length, i n is the number of such walks. C n is the set of 
n-step SSAWs whose last segment is at least one greater than the length of 
the third last segment and c n the number of such walks. From Figs. 2(a) - 
2(c), the section of the spiral from O to A belongs to and that from A to 
N belongs to C n _fc. The walks of the type shown in Figs. 2(a) - 2(b) can be 
uniquely decomposed in this way. Thus 

n— 1 n— 1 

Sn = ^Cn-k - ikin-k (6) 

fc=0 k=l 

with i = 1, c =0, Ci =1. The second sum corrects for a double-counting 
of walks of the type 2(c) by the first sum. A geometrical construction 



transforms walks of the type 2(c) into two connected walks from the set I. 
The walk in the set C n+ \ can be obtained by adding a step in a particular 
direction to each member of C„U I„. Thus 



-n+l 



C n + in (7) 



Now c n is equal to the number of partitions of the integer n. An example is 
given in Fig. (3). Fig. 3 (a) represents a particular partition of the number 12. 
A line at 45° divides the set of points into two classes, representing vertical 
and horizontal segments. There are two horizontal segments of length 2 and 
3 and two vertical segments of length 3 and 4. The resultant SSAW is shown 
in Fig.3(b). Each partition of the integer 12 can be represented by a pattern 
of points and each such pattern can be uniquely represented by a SSAW. All 
the spiral walks in the set C n can be enumerated in this way. 

Asymptotic expressions (n — > oo ) for the number of partitions of integers 



are known from the classic work of Hardy and Ramanujan||15||. One obtains 



e** 4 [1 + 4= + Of-)] (8) 



where c is a constant which can be determined. From (7), 



Substituting (8) and (9) into (6), one finds that the largest terms in the sum 
are in the vicinity of k = n/2. Expanding around the maximum and replacing 
the sum by an integral, one finally obtains 
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n~ exp[2 7r - ][l + — ] (10) 



Joyce [Pf] has obtained the complete asymptotic expansion for s n . Lin|T4" 
has obtained an expression for s n in the case of the triangular lattice. By 
comparing the expression (10) with that (Eq.(4)) valid for undirected and 
directed SAWs, one finds that the rotational constraint produces a new scal- 
ing form for the number of SAWs. Blote and Hilhorst |ll|] have also obtained 



an asymptotic form for the root mean square end-to-end distance of SSAWs. 
Again, the rotational constraint has a non-trivial effect on the scaling form. 



III. Spiral site lattice animals 

Two versions of the spiral LA have been studied : 'rooted', i.e., the origin 
from which the animal is grown is kept fixed and 'unrooted', i.e., the origin 
is not specified. In the latter case, two animals differing only in the choice 
of the origin are considered identical. The number of distinct rooted animals 
of a certain size is much greater than the number of unrooted animals of the 
same size. Bose and Ray|| have defined unrooted spiral LAs in d = 2 to 
be a subset of undirected and unrooted LAs with the proviso that there is a 
site of the cluster, namely, the origin to which every other site of the cluster 
is connected through at least one spiral path. In a spiral path, connection 
is either in the forward direction or in a specific rotational direction, say, 
clockwise. An example of a spiral site LA is shown in Fig.4. In dimension d 
I 2, a possible definition of the spiral constraint is: the projection of the path 
joining any site of the animal to the origin (the first site to be occupied), on 
any specified plane, say the xy plane, should have no anticlockwise turns. 



Other definitions of the spiral constraint have also been suggested pH[ . 

A variety of techniques can be used to study the configurational proper- 
ties of LAs. For spiral LAs, the method of exact enumeration of clusters upto 
a certain size has been used to obtain the values of A and the exponents 9 
and v in the scaling relations (4) and (5). Unlike in the case of SSAWs, the 
scaling relations (4) and (5) appear to give an adequate description of the 
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configurational properties of spiral LAs. The exact enumeration of animals 
of a certain size s is done on the computer using the well-known Martin's 
algorithm W7\. This algorithm generates undirected LAs and then each an- 



imal is checked for spiral connection. Knowledge of LA numbers of various 
sizes enables one to calculate the ratios g s +i/ g s for various values of s. From 
Eq.(4), in the large s limit, the plot of the ratio g s+ i / g s versus 1/s is that of 
a straight line. From the intercept and slope of the straight line, the growth 
parameter A and the exponent 9 can be determined. For spiral site LAs 
on the square lattice, exact enumeration data is available for size upto s = 
15|T2|. The values of A and 9 are determined as A = 3.002 ± 0.020 and 
9 = 0.69 ± 0.10. The radius-of-gyration exponent v = 0.50 ± 0.02, i.e., the 
LAs are nearly-compact. These values are different from the corresponding 
values in the cases of undirected and directed LAs. Thus spiral LAs belong 
to a new universality class. 

For the undirected and directed LAs, it is well-known that animals and 
trees (animals without loops) belong to the same universality class, i.e., loops 
have no effect on cluster statistics. This fact has been arrived at through 
detailed studies using techniques like both field-theoretic and position space 



RG, series-expansion and Monte Carlo enumeration [IS, [20], gTJ g^]. For spiral 
trees, the values of Ao, 9q and uq have been obtained as||18|| (from exact 
enumeration data upto s = 15) A = 2.123 ± 0.004, 9 = -1.318 ± 0.020 and 
uq = 0.67 ± 0.02. These values are different from the corresponding values 
for spiral LAs. Thus the spiral LAs and trees belong to different universality 
classes, i.e., loops have a non-trivial effect on spiral animal statistics. It has 
as yet not been possible to rigorously explain this fact but a special feature of 
spiral clusters is worth mentioning. For both ordinary and directed animals, 
an animal with loops becomes a tree as soon as a minimal number of bonds 
( for bond LAs) or sites (for site LAs) are removed but this is not always the 
case for a spiral animal. A tree obtained from a spiral animal with loops need 
not be a spiral tree. Also, in order to grow in the forward direction, a spiral 
animal should contain many loops. Without the loops, spiral connectivity 
is not possible in many cases. How these facts affect asymptotic cluster 
charateristics is still not understood. 

An interesting result associated with the LA problem is that of dimen- 



sional reduction. Parisi and Sourlasp3[ have shown that the asymptotic 



properties of LAs in d dimensions are the same as the critical properties of 
the Lee- Yang edge singularity (LYES) of the Ising model in d - 2 dimen- 
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sions. For the directed LA problem, the dimensional reduction is d - 1 . 
The LYES problem has to do with the fact that ferromagnetic Ising models, 
besides exhibiting singular behaviour at the critical point T = T c , in the 
presence of zero external magnetic field, show singular behaviour also at a 
critical value of imaginary magnetic field H = H (T) and for T i T c . The 
magnetiation has a branch point singularity of the form 

M — M ~ (H — HqY (11) 

where M is the magnetization at H . Consequently, the free energy F of the 
(d - 2)-dimensional Ising model has the same singular structure as that of 
the generating function of the d-dimensional undirected LA problem. From 
Eqs.(3) and (4), it follows that the LA generating function G(x) has a sin- 
gularity of the form (x — Xc) 6 " 1 where x c = 1/A. Since magnetization is 
the first derivative of the free energy F w.r.t. the magnetic field H, F has 
a power-law singularity of the type (11) with the exponent a + 1. From 
the equivalence of the free energy and the generating function in appropriate 
dimensions, one obtains 

0(d) = a(d- 2) + 2 (12) 

The Josephson scaling law (the hyperscaling relation) connects the singu- 
larity of the free energy F with the singularity of the correlation length £, 
(F ~ £,~ d )- One thus obtains 

v(d) = [<r(d-2)+]/(d-2) (13) 

Each undirected and unrooted LA of N sites gives rise to N distinct rooted 
site animals. Thus the animal number exponents, in the two cases, differ by 
1. Thus, from (12) and (13), one gets 

6 = (d-2)u (14) 

for rooted undirected animals. Similarly, one can show that 

= (d-l)u ± (15) 

for rooted directed animals. The exponent u± is associated with the radius 
of gyration of the directed animals perpendicular to the preferred direction 
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of growth. We have seen that the LA exponents 9, v are related to the 
LYES exponent a. The value of a can be calculated exactly for d = 
(point particle) and d = l.Thus from Eqs.(12) and (13), the exponents 9, v 
can also be determined exactly in appropriate higher dimensions (upto d = 
3). The directed LA problem is exactly solvable in both 2d and 3df2~5|.So 
exploiting the fact that directed animal exponents in d dimensions are related 
to the ordinary animal exponents in (d + 1) dimensions and to the LYES 
exponents in (d-1) dimensions, 8, v are exactly known for ordinary animals 
in d = 4 and a, the LY exponent, in d = 2. The results for a(d), 6(d), 
v(d) are: a(0) = — l,cr(l) = —1/2, cr(2) = —1/6, for ordinary animals 
6(2) = 1,6(3) = 3/2, 6(A) = 11/6, u(3) = 1/2, v(A) = 5/12 and for 
directed animals 6(2) = 1/2, 6(3) = 5/6, v±(2) = 1/2, z/j_(3) = 5/12. 

To study dimensional reduction in the spiral animal problem, Bose et 
al|]rm have considered only spiral trees. Both exact enumeration and MC 
simulation have been performed. The spiral trees considered are rooted and 
the rotational constraint is operative only in the xy plane. The radius of 
gyration has two components R s , p i in the xy plane and R s ,± in the direction 
perpendicular to the xy plane. The numerical data obtained are consistent 
with the following results: 

#o = (d-A)u 0yPl (d = 2) 

6 = (d-A)v , ± (d>2) (16) 

Eq.(16) shows that a dimensional reduction by 4 occurs in the spiral animal 
problem. The dimensional reduction by the factor two, in the case of ordi- 
nary animals, occurs due to a hidden supersymmetry [p3|, [2B|. The first few 



applications of supersymmetry in condensed matter and statistical physics 
include the following problems: the random field model, the linear polymers 
(SAWs on a lattice), the branched polymers (LAs) and electron localization. 
These models can be suitably recast into forms which manifest supersymme- 
try. Introduction of an effective supersymmetry simplifies technically com- 
plicated calculations. Also, as in the cases of branched polymers and random 
field models, which are closely related, the supersymmetry can be shown 
to lead to an effective reduction in the number of degrees of freedom. For 
the directed animal problem, the dimensional reduction by 1 occurs because 
the disorder or cellular automaton condition is satisfied. Field theories have 
been formulated for both ordinary and directed animal problems, specially, 
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with the view to understand the mechanism of dimensional reduction. For 
spiral LAs, a complete theory is yet to be worked out. A simple example of 
reduction of effective dimensionality due to rotational constraint, is provided 
by the quantum mechanical motion of a charged particle in a uniform mag- 
netic field. If the field is perpendicular to the xy plane, then the degrees of 
freedom in this plane get quantized. The resulting discrete spectrum has no 
contribution to the density of states of low energy excitations. This leads to 
a reduction by 2 in the effective dimensionality. One can argue, in analogy, 
that for the spiral animal problem also, the rotational constraint is respon- 
sible for a dimensional reduction by 2. A further reduction by 2 occurs due 
to the conventional Parisi-Sourlas mechanism. 

Recently, the problem of rooted spiral site trees embedded in a triangular 



lattice has been studied by exact enumeration||27|| . It will be of interest to 
verify whether the dimensional reduction by 4 occurs in this case also. The 
universality of the exponents 9 and v imply that they have the same values 
as long as the dimension d of the underlying lattice is the same. Also, the 
values remain unchanged if one considers bond LAs instead of site LAs. These 
facts are well-established for both the undirected and directed LA problems. 
Preliminary results in the case of the spiral LA problem, however, indicate 
that the values of 9 and v are different for the square and triangular lattices. 
If the result is proved to be true, then one has an example of the breakdown 
of universality in the case of the LA problem. 

We next consider the configurational statistics of spiral site LAs with c (c 
= 0,1,2,...) number of loops. The corresponding growth parameter and the 
exponents are A c , 9 C and v c respectively. The number of loops c in an animal 
of s sites has been determined through the relations: 

2 

6 = 1/2 53 iq(i) ,c = b - s + 1 (17) 

i=l 

where b is the number of bonds in the animal, z is the coordination number of 
the lattice and q(i) is the number of sites in the animal which have i nearest- 
neighbour occupied sites. In the case of undirected and directed LAs, A c and 
u c are independent of c and equal to A and u , the respective values for trees. 
The animal number exponent 9 C is related to that for trees (9q) through the 
relation [2|, |2| 

9 C = 9 - c (18) 
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A rigorous proof of (18) has been given for undirected weakly-embeddable 
animals ( bond animals) in the square lattice [Q. Universality demands 
that the same relation be true for strongly-embeddable (site) animals. For 
spiral LAs, no rigorous result is available but calculations based on exact 
enumeration data suggest the relation 

e c = 9 - ac (19) 



where a = 1.5 ± 0.2.Whittington et al |28[ have derived the inequality 

0c > 0c+i > C - 1 (20) 
The inequality is derived from another inequality involving LA numbers 

9s-3,c < 9s, c+i < 2dsg S:C (21) 
where g s c ~ A* s~ 6c . For spiral site animals, an inequality similar to the first 



half of the inequality in (21) can be proved |T8[. Let the vertices or sites of the 
spiral animals have coordinates (xi,yi), i=l,2,....,n.Following Whittington et 
al |pS|1 , the top (bottom) vertex is defined as the vertex having maximum 
(minimum) x coordinate and, in case of ambiguity, the vertex in this sub- 
set having maximum (minimum) y coordinate. Consider any c-loop animal 
belonging to the set W SjC of animals with s sites and c loops. The number 
of animals is g SjC . Let t be the top vertex of the animal with coordinates 
(x t , yt)- The top vertex can be approached in three ways: (i) from the vertex 
(xt — 1, yt) if it is occupied, (ii) from the vertex (x t , y t — 1) if it is occupied, 
(iii) from both the vertices (x t — 1 , y t ) and {x t , y t — l) if both are occupied. 

In cases (ii) and (iii), add four vertices vi, v%, V3, V4 with coordinates (x t + 
l,y t ), {x t + 1, y t - 1), (x t + 2, y t ) and (x t + 3, y t ) respectively. For lattice 
site animals, addition of a new site implies that the bonds between this site 
and its nearest occupied neighbours automatically exist. Addition of four 
new sites leads to an increase in the number of bonds by five. The number 
of loops therefore increases by one. The new animal obtained is also spirally 
connected and so belongs to the set W s+ ^ c+ i of animals with c+1 loops. In 
case (i), there are two possibilities, either the vertex (x t , yt — 1) is occupied 
or unoccupied. If the first possibility is true, add four vertices with the same 
coordinates as in the cases of (ii) and (iii) to get a spiral animal of s + 4 sites 
and c+1 loops. If the second possibility is true, add vertices with coordinates 
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(x t + l,yt), (x t + 2, y t ), (x t + 1 , y% - 1) and (x t + 2, y t - 1) generating a 
spiral animal of s + 4 sites and c + 1 loops. Thus for each animal belonging 
to the set W s>c of animals, an animal belonging to the set W s+4jC+ i can be 
generated ,i.e., g SjC < g s +4, c +i or, g s _ AyC < g S)C+ i. 

The proof given by Whittington et al[28| for the second half of the in- 
equality fails in the case of spiral animals. For undirected weakly embeddable 
animals, an upper limit to g s>c is 2dsg s<c ^i as there are s vertices and at each 
vertex an edge can be added in 2d ways, addition of an edge (bond) im- 
plying an increase in the number of loops by 1. Thus g S}C < 2dsg S}C -i, or, 
9s,c+i < 2dsg SjC , the second half of the inequality in (21), is obtained. In the 
case of spiral lattice animals, animals with c loops can be generated not only 
from spiral animals with c-1 loops but also from undirected animals with c - 
1 loops, which on the addition of an edge( bond) become spirally connected 
with c loops. Thus a simple relation like g SjC < 2dsg s ^ c _i cannot be written 
down for spiral lattice animals. 

For spiral LAs, no rigorous results are available. This is in contrast to 
the case of directed LAs. In this context, we mention a directed compact 
site LA problem on the d-dimensional hypercubic lattice |3T|. It has been 
rigorously established that the model is equivalent to the (i) infinite-state 
Potts model and (ii) the enumeration of (d-l)-dimensional restricted parti- 
tions of an integer. Thus, the directed compact LA problem can be solved 
exactly in d = 2,3 using known solutions of the enumeration problem men- 
tioned. Also, knowledge of the infinite-state Potts model solution leads to a 
conjectured limiting form, for d ^3, of the generating function of restricted 
partitions. The latter problem is still unsolved in number theory. The ex- 
amples of SSAW and the directed compact site LA problem show that the 
well-known mathematical theory of partition of integers can be used to derive 
rigorous results for lattice-statistical models. 



IV. Spiral site percolation 

Santra et al[^, |33| have formulated a new type of percolation model known 
as spiral percolation in which a rotational constraint is operative due to 
which each step in a percolation path is either in the forward direction or 
in a specific rotational direction, say, clockwise. Consider spiral site per- 
colation on a lattice in 2d of size LxL. The spiral percolation threshold is 
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determined by the binary search method One starts from a central site 
of the lattice called the origin. Initially, from the origin, one can proceed in 
any one of the four possible directions. The nearest-neighbour (n.n.) sites in 
these directions are occupied with probability p by using a random number 
generator. With each occupied site a variable JVISIT(IV,INOW) is asso- 
ciated. INOW is the site index and IV the direction index. The sites of 
the lattice are numbered according to a particular sequence; the site index 
is the number associated with a site. A site can be reached from four di- 
rections,south,west, north and east. The corresponding direction indices are 
1,2,3 and 4 (Figure 5(a)). The rotational constraint implies the following: 
if a site has direction index 1, then occupy with probability p, only those 
of the n.n.s which are towards the north and east of the site. Initially, the 
JVISIT variable is assigned the value zero for all sites. As soon as a site is 
occupied, the corresponding JVISIT variable is given the value IV where IV 
is the direction index corresponding to the direction from which the site is 
occupied. All growing sites are put on a list and the walk-search is carried 
out for each. The growth of a cluster of occupied sites stops only when all the 
perimeter sites are unavailable for occupation. Figure 5(b) shows a typical 
cluster grown obeying the rotational constraint. Due to the nature of the 
constraint, loops are an essential feature of the growing spirally connected 
cluster. In the computer algorithm, the possibility of loop formation while 
growing a cluster is taken care of in the following manner. Refer to Figure 
5(b). The site i is the origin. In the first step when j is occupied, JVISIT(iJ) 
=1. Next time j is approached (from site m,say), the computer algorithm 
checks whether JVISIT(4,j) is zero. If it is non-zero then the site j has been 
approached from the east previously and so cannot be occupied again. For 
the particular cluster in Figure 5(b), since JVISIT(4,j) is zero, the site j is 
reoccupied and JVISIT(4,j) assigned the value 4. Thus loop formation and 
continuation of cluster growth are possible. 

Following the binary search method, one starts growing a cluster with 
a particular value of p,say p . If the cluster grown spans (does not span) 
the lattice in either the east-west or north-south direction, p is decreased 
(increased) by a small amount. The same random number sequence is then 
used to get estimates Pi(L) and P2^L) which bound an interval containing 
the true threshold value p c {L). By successive binary chopping of this interval 
one can determine p c (L) with a specified accuracy. The whole process is then 
repeated N times (N x L x L is of the order (10 6 — 10 7 )) using different 
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random number sequences. The average value < p c (L) > of all the estimates 
obtained is taken as an estimate for the percolation threshold. 

We now describe the scaling behaviour of the percolation clusters in the 
vicinity of the percolation threshold. In the case of undirected percolation, 
the most widely used method of calculating the scaling behaviour is to fill up 
a large lattice by clusters of various sizes s with site occupation probability 
p||. The cluster size distribution is given by n s (p), the number of clusters of 
size s per site of the lattice. The probabilty that a site belongs to a cluster of 
size s is sn s (p). The probability that an arbitrary site belongs to any cluster 
is equal to the probability that it is occupied. Thus 

J2 sn s (p) = p (p < p c ) (22) 

s 

Y sn s (p) + Poo(p) = p{p > p c ) (23) 

s 

The sums in (22) and (23) run over all possible finite cluster sizes. Poo{p) is 
the probability that an arbitrary site (occupied or not) belongs to the infinite 
cluster. The percolation probability P(p) is defined as the probability that 
an occupied site belongs to the infinite cluster, i.e., 

P{p) = P 00 / p = 1 -J2 sn s (p) (24) 

s 

Let W = (sn s )/ (J2 S sn s) be the probability that the cluster to which an 
arbitrary occupied site belongs contains exacly s sites. The average cluster 
size S av is then given by 

Sav =Y sW s = H (SV)/(E Sn s) (25) 
s s 

The linear size of the finite clusters, below and above p c , is characterised 
by the correlation length £. The correlation length is defined as the mean 
distance between two sites on the same finite cluster. As the site occupation 
probability p reaches the percolation threshold p c , the various moments of 
the cluster size distribution n s (p) given by 

£s fe n s (p), k = 0,1,2,... (26) 
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(the prime denotes that the largest cluster is excluded from the sum) be- 
come singular with characteristic critical exponents. The critical exponents 
a, j3, 7, 5, v are defined through the relations 

k = : (£ n s ) sing ~ | p - p c I 2 "" (27) 

s 

I 

k = 1 : P(p) oc C£sn s ) sing ~ (p-Pc)" (28) 

s 

/ 

k = 2 : S av oc ^ s 2 n s ~ | p - p c (29) 

s 

e=|p-p c |- / (30) 

Y,sn s ( Pc )(l -e- hs ) ~ fcV« (31) 

The subscript 'sing' refers to the leading non-analytic contribution. The 
above equations are valid in the limit in which the fictitious field h and the 
quantity | p — p c \ both go to zero. 

We now consider the case of spiral site percolation. The scaling behaviour 
of clusters in the critical region is studied following the single-cluster growth 
method in which clusters are grown singly starting from a fixed origin instead 
of filling up the whole lattice by clusters of various sizes. The computer 
algorithm has already been described. The cluster size distribution is now 
defined as 

P s (p) = Ns/Ntot (32) 

where N s is the number of clusters of size s in a total number N tot of clusters 
generated. The various moments of P s {p), J2' s $ k Ps{p) becomes singular as 
p — > p c . The prime in the summation symbol indicates that the largest 
cluster is to be left out of the sum. The average cluster size x is given by 

X~E^b) (33) 

s 

which diverges as x ~ I V ~ Pc |~ 7 as p — > p c . In the single cluster growth 
method, the cluster configurations are rooted at the origin and the first, not 
the second, moment of the cluster size distribution defines the exponent 7. In 
fact, any exponent which corresponds to the kth moment of the cluster size 
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distribution in the case of the lattice-filling method, is obtained from the (k- 
l)th moment, one order less, in the case of the single cluster growth method. 
In the case of directed and spiral percolation, where there is a directionality 
constraint on the percolation process, the single cluster growth method is the 
appropriate one for the study of critical behaviour. 

The percolation transition is a second-order phase transition and exhibits 
the critical phenomena characteristic of second-order thermodynamic phase 
transitions like liquid-gas and paramagnet-ferromagnet transitions. The ma- 
jor characteristic feature, which we have already discussed, is that different 
cluster-related quantities exhibit power-law behaviour as p — + p c . One of 
these, the percolation probability P(p) (Eq.(24)) serves as the order parame- 
ter of the transition. The order parameter has a non-zero value in one phase 
(p I Pc) and is zero in the other phase (p j p c ). As in the case of ther- 
modynamic phase transitions, the critical behaviour is universal, i.e., the 
magnitude of the critical exponents depends only on the dimensionality of 
the lattice and is the same for both site and bond percolation. The critical 
exponents are not independent (only two are so) and one can write down 
'scaling' relations connecting different critical exponents. The scaling rela- 
tions can be determined on the basis of the scaling hypothesis according to 
which in the vicinity of p = p c , the cluster-related quantities become gener- 
alised homogeneous functions (GHF). A function F(x,y) of two variables is 
a GHF if it satisfies the property 

F(x,y) = x A f(y/x B ) (34) 

where the function f is the appropriate scaling function and is the function of 
a single variable. The form (34) usually holds true when both the variables x 
and y tend to zero. The definition (34) of a GHF can be generalised to more 
than two variables. A second-order thermodynamic phase transition has the 
feature that fluctuations occur at all length scales, i.e., the correlation length 
goes to infinity. In the case of percolation transition also, the correlation 
length goes to infinity. The infinite percolation cluster is a self-similar object 
with fractal dimension less than the Euclidean dimension of the space in 
which the cluster is embedded. 

We now describe the approximate methods for calculating the critical ex- 
ponents and determining the scaling behaviour of the cluster-size disribution 
function. The methods normally used are Monte Carlo (MC) simulation, 
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finite-size scaling, series-expansion and renormalization group. We consider 
the case of spiral site percolation on the square and triangular lattices and 
describe the results obtained by Santra and Bose[B2], B3. In MC simulation, 



a particular lattice size is considered. The value of the percolation threshold 
p c is determined by the binary search method. Ten thousand cluster are 
generated by the single cluster growth method. The growth of a cluster is 
stopped when there is no further site available for occupation. The quanti- 
ties, the average cluster size x (Eq.33 ), the second moment x' °f the cluster 
size distribution P s (p) and the correlation length £ are determined as a func- 
tion of | p — p c | . As the site occupation probability p tends towards the 
percolation threshold, \i x' an d £ diverge with the exponents 7, 7' and v. 
We first quote the results for the square lattice. The value of p c is p c = 0.711 
± 0.001. The critical exponent 7 is obtained from the slope of logx versus 
log I p — p c I and is given by 7 = 2.19 ± 0.07.The exponents 7' and v are 7' 
= 4.51 ± 0.16 and v — 1.16 ± 0.01. The exponent v can also be determined 
from the scaling relation v = (7' — 7) /2 as v = 1.18 ± 0.23. This value is in 
agreement with the value obtained through direct measurement. For the tri- 
angular lattice, p c = 0.667 ± 0.001, 7 = 2.079 ± 0.062, 7' = 4.376 ± 0.199 
and v = 1.012 ± 0.025. 

True critical behaviour occurs only in the limit of infinitely large lattices. 
An estimate of the critical exponents can, however, be obtained from the 
studies of finite systems by assumung the finite-size scaling hypothesis [§] 
which leads to the formula 

A = L-*/ v F\(p-pc)L l /»] (35) 

where A is a quantity which becomes critical, A ~ | p — p c \ x as p — > p c - The 
function F is a suitable scaling function. At p = p c , the quantity A varies 
as L~ x l v . This result can be used to determine the exponents 7/z/, 7'/V by 
calculating the first and second moments x an d x' respectively, of the cluster 
size distribution. The percolation threshold p c = jp c (L)^, the average being 
taken over a large number of estimates of p c for a L x L lattice. The largest 
cluster spanning the lattice is of size Sl and is a fractal with fractal dimension 
df defined by 

S L ~ L*i (36) 
The fractal dimension df can again be written as 

d f = d - (3/u (37) 
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where (3 is the exponent associated with the order parameter. From (36), 
the fractal dimension df can be determined and from (37) putting d = 2, the 
value of (3/v can be calculated. The measurement of various quantities is 
done for lattices of various sizes. The slope of the straight line logS^ versus 
logL gives d f = 1.965 ± 0.009 and so (3/v = 0.043 ± 0.009. From (35), the 
plot of logA versus logL at p = p c gives the exponent x/v. Thus for A = x, 
the average cluster size, the exponent j/v — 2.01 ± 0.06. For A = x', the 
second moment of the cluster size distribution, the exponent 7'/ v = 4.05 ± 
0.13. For the triangular lattice, the various quantities have been determined 
as d f = 1.965 ± 0.008, (3/v = 0.035 ± 0.008, 7/1/ = 1.867 ± 0.028 and 
i/v = 3.829 ± 0.062. 

We now describe the method of series expansion. The probability p that 
the origin from which clusters grow is occupied can be written as a sum over 
all finite clusters that start from it (for p j p c ), i.e., 

P = T,9stP s (l-pY = Y,P SD M, q = 1-P (38) 

s,t s 

where g s t is the number of clusters or LAs of s sites and t perimeter sites and 
-D s (g)'s are the perimeter polynomials. The animals are rooted at the origin. 
The average cluster size x is the first moment of the cluster size distribution. 

X = sg st p s (l-pY Hp - Pc I" 7 (39) 

s,t 

as p — > p c . The clusters rooted at the origin obey the same rotational con- 
straint as in the percolation process. The spiral LAs, as already mentioned 
in Section III, have been extensively studied by Bose et al || [H| |18| . The 
major focus in these studies has been on spiral site LAs with no fixed origin. 
If the origin is not kept fixed, the number of perimeter sites is dependent on 
the position of the origin. In spiral percolation, the percolation clusters are 
rooted at the origin so the spiral LAs in the series-expansion method have to 
be rooted, i.e., have a fixed origin. Now each distinct LA configuration has a 
unique value for the number of perimeter sites. The rooted spiral site LAs on 
the square lattice have been enumerated exactly by Santra and Bose[[33| for 
size upto 13. Thus, series expansion for the average cluster size x ( ec L-(39)) is 
exactly known, when expanded in powers of p, upto the term containing the 
power p 13 . A standard trick [|33[ enables one to extend the series by one more 
power. Thus the coefficients in the expansion of x — 1 + J2 r b r P r (a common 
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factor of p has been taken out), are exactly known for r = 1 to 13. As p — > p, 
X has the power-law behaviour x ~ ijPc — p)^ = Pc 1 (1 ~~ v/Vc)^ 1 ■ On bino- 
mial expansion, the ratio of coefficients b^+i/bk = l/p c (l + (7 — l)/(k+ 1)). 
The plot bk+i/bk versus l/(k+l) should be a straight line. The intercept 
of the line on the y-axis gives l/p c and from the slope, knowing p c , the 
exponent 7 can be determined. A more sophisticated method of obtaining 
the value of 7 is that of Pade approximants. Using this method, the value 
of 7 for spiral site percolation on the square lattice has been obtained as 
7 = 2.167 ± 0.004. For lack of a sufficiently long series, the series-expansion 
method has not been applied to the triangular lattice. The series-expansion 
method can be applied to generate the series for other percolation-related 
quantities and thereby obtain the values of the associated exponents. 

We now discuss the scaling hypothesis for percolation clusters. From 
Eq.(34), one finds that if the function G(z) = F(x,y)/ x A is plotted as a 
function of z = y/x B , all the data (for different values of x) fall on a single 
curve described by the function f(z). On the other hand, if Eq.(34) were not 
true, the data, for different values of x, fall on separate curves. The collapse 
of data on a single curve serves as a verification of the scaling hypothesis. 
The cluster size distribution function for sufficiently large s has the scaling 
form in the critical region given byfl33 



Ps(p) = s~ T+1 fi(p- Pc )s°] (40) 

where both p - p c and 1/s go to zero as p — > p c . The exponents r and a are 
related to the exponents 7 and (3 through the relations 

7 = (3-r)/<7 ,(3 = (t-2)/<t (41) 

A verification of the scaling function form (40) is possible by plotting P s (p) / P s (p c ) 
against (p— p c )s°". If the scaling form is true then, for sufficiently large clusters 
and for different values of p, the data should collapse onto a single curve. For 
the square and triangular lattices, the scaling form (Eq.(40)) for spiral site 
percolation has been verified with a = 0.446 (square lattice) and o = 0.473 
(triangular lattice). 

The values of the critical exponents obtained for spiral site percolation in 
2d are different from the values of critical exponents in the case of undirected 
percolation. Spiral percolation is similar in nature to directed percolation, 
since for both the models a directional constraint is operative. In the first 
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case, the directional constraint is rotational in nature and in the second case 



percolation occurs only in certain specific directions. Dhar and Barmap4 
have studied directed site percolation on the square lattice using MC simu- 
lation. A comparison with their results shows that the average cluster size 
exponent 7 = 2.19 ±0.03 (directed percolation) is the same as that obtained 
in the case of spiral percolation, 7 = 2.19 ± 0.07, using MC simulation. 
Thus the average cluster size diverges with the same exponent irrespective 
of the nature of the external constraint. One crucial difference between spiral 
and directed percolation is that in the first case, percolation clusters grow 
isotropically whereas in the second case, the clusters grown are anisotropic 
in nature. For the first case, there is thus only one correlation length which 
diverges with the exponent v as p — > p c , whereas, in the second case, there 
are two correlation lengths, one parallel and the other perpendicular to the 
preferred direction, which diverge with the exponents un and u±_ respectively 
as p -> p c . 

One important quantity in the percolation problem is the cluster exter- 
nal perimeter or 'hull' of the large clusters. The hull of a cluster is the 
continuous path of occupied sites at the external boundary of the cluster. 
The hull has a significance of its own apart from being a part of the per- 



colation cluster. Sapoval et al|35[ have shown that the diffusion front arising 
out of diffusion from a source has a fractal structure that is related to the 
hull of percolation clusters. The percolation cluster hull in the case of undi- 
rected percolation exhibits scaling behaviour characterised by critical expo- 
nents which have values different from those of analogous quantities in the 



case of percolation clusters. Santra and Bose |36|] have studied the scaling 
behaviour of the percolation hull of large clusters in the case of spiral site 
percolation on the square and triangular lattices. The scaling behaviour is 
found to be different from that of hull in the case of undirected percolation. 
In the latter case, Saleur and Duplantier[|37| have proved that the fractal 
dimension, dn, of the hull is given by 

4=1+ l/v (42) 

where v is the correlation length exponent. Eq.(42) is an exact result. A 
similar relation cannot be proved in the case of spiral site percolation. 
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V. Concluding Remarks 



The rotational constraint, we have seen, has a non-trivial effect on the statis- 
tics of lattice models like SAW, LA and percolation. The undirected SAW 
and LA problems have been extensively studied and some fairly rigorous 
results are available. For undirected percolation in 2d, the values of p c are 
known exactly in some cases. The critical exponents are, through some exact 
mappings, conjectured to be known exactly 0]. For the directed SAW, LA 
and percolation problems, many exact /analytical results are available. In 
the presence of a rotational constraint, however, the only rigorous results 
that are known are for the SSAW model. The spiral LA and spiral perco- 
lation problems have been studied through numerical methods. Derivation 
of numerical results for the models is a challenge for both mathematicians 
and physicists. The computational methods in the case of spiral models are 
also more difficult to implement than in the case of undirected and directed 
models. 

Lattice models have extensive applications in all branches of science. The 
SAW model describes a correlated walk, i.e., a walk with memory. The 
knowledge gained from the study of the model has been found to be useful 



in problems like correlations in nucleotide sequences ||39|| and the folding of 



proteins [[EJ. Several examples of the various applications of the lattice mod- 



els can be found in the Refs.f 1-5]. Examples of the spiral structures and the 
rotational constraint can be found in Ref. [41]. A few examples are: spiral 
galaxies, the pattern of shoot arrangement in plants, scales on a pineapple 
or on a pine cone, the Belousov-Zhabotinsky reaction, plasma in a stochastic 
magnetic field, polymers with spiral structure, bacterial colonies which ex- 
hibit chiral morphology ||42|| (the colony consists of twisted branches all with 
the same handedness) etc. It is to be hoped that further studies will be un- 
dertaken to understand the effect of rotational constraint on the properties 
of model systems. 
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Figure Captions 



Fig.l A spiral self- avoiding walk with an outward spiralling part and an 
inward spiralling part; 'O' denotes the origin. 

Fig. 2 Three different spiral self- avoiding walks with origin O and end-point 
N. The broken line divides the walk into two parts. 

Fig. 3 (a) A partition of the number 12 into the integers 4,4,2 and 2. (b) the 
spiral self- avoiding walk corresponding to the set of points in (a). 

Fig. 4 A spiral lattice animal of 10 sites on a square lattice grown from the 
origin X. 

Fig. 5 (a) Direction indices of a site corresponding to the four different di- 
rections from which the site can be reached, (b) an example of a spiral 
site cluster on the square lattice. The arrows on the bonds indicate the 
allowed spiral directions of flow from site i. 
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